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Abstract 

We have compared the 14-year record of satellite de- 
rived tropical tropospheric ozone columns (TTOC) from the 
NIMBUS-7 Total Ozone Mapping Spectrometer (TOMS) to 
TTOC calculated by a chemistry-transport model (CTM). An 
objective measure of error, based on the zonal distribution 
of TTOC in the tropics, is applied to perform this compari- 
son systematically. In addition, the sensitivity of the model 
to several key processes in the tropics is quantified to se- 
lect directions for future improvements. The comparisons 
indicate a widespread, systematic (20%) discrepancy over 
the tropical Atlandc Ocean, which maximizes during austral 
Spring. Although independent evidence from ozonesondes 
shows that some of the disagreement is due to satellite over- 
estimate of TTOC, the Atlantic mismatch is largely due to 
a misrepresentation of seasonally recurring processes in the 
model. Only minor differences between the model and ob- 
servations over the Pacific occur, mostly due to interannual 
variability not captured by the model. Although chemical 
processes determine the TTOC extent, dynamical processes 
dominate the TTOC distribution, as the use of actual mete- 
orology pertaining to the year of observations always leads 
to a better agreement with TTOC observations than using a 
random year or a climatology. The modeled TTOC is re- 
markably insensitive to many model parameters due to effi- 
cient feedbacks in the ozone budget. Nevertheless, the simu- 
lations would profit from an improved biomass burning cal- 
endar, as well as from an increase in NO x abundances in free 
tropospheric biomass burning plumes. The model showed 
the largest response to lightning NO x emissions, but system- 
atic improvements could not be found. The use of multi-year 
satellite derived tropospheric data to systematically test and 
improve a CTM is a promising new addition to existing meth- 
ods of model validation, and is a first step to integrating tro- 
pospheric satellite observations into global ozone modeling 


Correspondence to: W. Peters (W.Peters@phys.uu.nl) 


studies. Conversely, the CTM may suggest improvements to 
evolving satellite retrievals for tropospheric ozone. 
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1 Introduction 

The availability of trace gas measurements in the tropics has 
increased considerably by a large number of intensive cam- 
paigns (LBA, PEM, TRACE, INDOEX, SAFARI) and the 
expansion of the number of measurement stations in the trop- 
ics. The interest in the tropics is based on the active pho- 
tochemistry, rapid convective mixing along with lightning 
discharge, high humidity, widespread biomass burning, and 
growing industrial emissions. All these processes influence 
the concentration of oxidants such as ozone and the hydroxyl 
radical. The latter plays a crucial role in the removal of bio- 
genic emissions and anthropogenic pollution from the atmo- 
sphere. 

Observations from satellites play an increasingly impor- 
tant role in monitoring of the chemical composition of 
the troposphere. The Total Ozone Mapping Spectrometer 
(TOMS) instrument aboard the NIMBUS-7 and Earth Probe 
(EP) Satellites (McPeters and Labow, 1996) has been used 
to retrieve tropospheric column ozone (TTOC) (Fishman and 
Larsen, 1987; Kim et al., 1996; Ziemke et al., 1998; Hudson 
and Thompson, 1998), while the GOME instrument aboard 
ERS-2 has been used to retrieve columns of tropospheric 
NO 2 , SO 2 (Eisinger and Burrows, 1998; Burrows et al., 
1999), and BrO (Richter et al., 1998). In the near future, 
SCIAMACHY (Burrows and Chance, 1992) will expand this 
suite of trace gases with N 2 O, CO 2 , CO and CH 4 . These last 
two trace gases are also the focal point of the MOPITT in- 
strument (Pan et al., 1998) aboard the EOS-TERR A satellite, 
launched in 2000. The most important advantages of satel- 
lite observations are the global coverage and the continuous 
monitoring on time scales from days to years, and sometimes 
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even decades. For historically data sparse regions such as the 
tropics, tills will contribute significantly to our understanding 
of atmospheric processes and help detect changes. 

In addition, satellite measurements can easily be combined 
with chemistry-transport models (CTM) and global circula- 
tion models (GCM), since spatial and temporal scales are of- 
ten comparable. There are several ways in which these prop- 
erties can be exploited in a quantitative way. This includes 
for instance assimilation of satellite measurements into mod- 
els to constrain chemistry calculations (e.g. .Teuken et aL, 
1999) and provide boundary conditions. Also, satellite mea- 
surements can be used as additional constraints in inversion 
problems, where the global distribution of a tracer is used to 
infer the distribution and magnitude of its sources and sinks 
(e.g. Houweling et al., 1999). Vice versa, model calculations 
can be used in satellite retrieval algorithms to provide a ’first 
guess’ field or to provide missing parameters needed in the 
retrieval (e.g. Chance et al., 2000). The growing amount of 
satellite observations and expansion of remotely sensed trace 
species will intensify the use of these observations in model- 
ing studies. 

From the several tropospheric satellite observations avail- 
able at the moment, the TTOC record derived from the 
NIMBUS-7 observations is the most extensive. It spans the 
complete NIMBUS-7 observational period (1979-1992) on 
a high spatial (2x1 degrees) and temporal (twice monthly) 
resolution. This has enabled analysis of seasonal cycles 
(Martin et al., 2000), trends (Thompson and Hudson, 1999) 
and interannual variability (Ziemke et al., 1999) in tropical 
ozone. Moreover, this dataset has shown to be very useful 
in studying both local scale features such as the effect of El 
Nino-Southern Oscillation (ENSO) related fires over Indone- 
sia (Chandra et al, 1998; Thompson et al., 2001), as well as 
large scale phenomena such as the evolution of the South At- 
lantic ozone maximum (Thompson and Hudson, 1999) and 
the changing Walker circulation during ENSO (Ziemke et al., 
1999). Neverthless, most modeling studies have only used 
the TTOC data for a qualitative comparison, while ozone 
sondes (Logan, 1999)) were used for a quantitative compar- 
ison. Here, we explore how satellite derived TTOC can be 
used to validate and improve our CTM. 

Previously, we have tested our model against in situ data 
from sondes and ground-based stations (Lelieveld and Den- 
tener, 2000; Marufu et al., 2000; Peters et al., 2001). Com- 
parison to the climatology of Logan (1999) revealed an ap- 
parent problem in the representation of the zonal ozone dis- 
tribution in the tropics, which shows a maximum of TTOC 
over the Atlantic Ocean and a minimum over the Pacific 
Ocean. This distribution is often called the ’zonal wave one 
pattern’ (Kim et al., 1996), and is a typical example of a large 
scale feature in tropical ozone that can be well studied with 
a combination of global models and satellite observations. 
Therefore, we present a quantitative comparison with TTOC 
derived with the Modified Residual (MR) algorithm of Hud- 
son and Thompson (1998). We have developed a measure 
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of error to consistently perform this comparison. This mea- 
sure of error allows us to recognize systematic shortcomings 
in our model on seasonal and interannual time scales, and 
study the influence of meteorology and photochemistry on 
our model results. Moreover, it helps to identify processes 
that should be considered for improvement and allows us to 
quantify the sensitivity of our model results to several pro- 
cesses such as biomass burning and lightning NO x emis- 
sions. It is important to stress that these techniques are not 
used to ’tune’ the model to observations. The processes that 
are most likely to cause the observed differences should be 
carefully studied in future work, preferably in a three dimen- 
sional perspective using ozone sondes. Our work presented 
in Section 8 should be seen as an exploration of model sen- 
sitivities providing guidance to actual improvements. 

We start in Section 2 with an introduction of the CTM that 
was used in this study, and a discussion of the photochem- 
istry, biomass burning and lightning NO x emissions. Next, 
we give some basic information on the TTOC data used (Sec- 
tion 3). In order to compare model and measurements, we 
introduce a measure of error in Section 4, which will be the 
basis for the results in this paper. We apply this measure of 
error to firstly identify the magnitude, location, and season- 
ality of the discrepancy between model and data (Section) 5, 
after which we will look at the roles of meteorology (Section 
6) and photochemistry (Section 7) in our model. The latter 
is done by looking at the spatial and temporal distribution of 
photochemically marked ozone tracers. The processes that 
are identified as most likely candidates to explain the dis- 
crepancy are further studied in Section 8, through a series of 
sensitivity experiments in which changes in emission mag- 
nitude and distributions are applied. The results from these 
experiments are illustrated with several ozone sondes in Sec- 
tion 9. Finally, we discuss our results in Section 10, and give 
a short summary in Section 11. 

2 Model Description 

We use the chemical tracer transport model version 3 (TM3) 
to simulate transport and photochemistry of 48 trace gases 
within a global domain. The TM3 model has a horizontal 
resolution of 5° longitude x 3.75° latitude and uses 19 lay- 
ers with a hybrid of sigma-pressure coordinates in the verti- 
cal, with the model top at 10 hPa. Approximately five layers 
constitute the boundary layer, nine the free troposphere, and 
five the lower and middle stratosphere. Transport in the off- 
line model is prescribed using ECMWF data (Gibson et al., 
1997). 

The model uses six-hourly, 3-dimensional fields of geopo- 
tential height, horizontal mass-fluxes, temperature, specific 
humidity and cloud information (liquid water and ice wa- 
ter content), as well as 2D information on cloud cover and 
precipitation (large-scale and convective). The subgrid scale 
transport by turbulence and convection is calculated from 
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Table 1 . Processes included in the model and references to the original work 

process Based on Ref-* 

Russel and Lemer (1981) 


Advection 

Convection (subgrid) 
Turbulent mixing (subgrid) 
Stratospheric boundary 
Photolysis 

Dry deposition 
Wet deposition 
Biomass burning emissions* 
Lightning NO x emissions* 

Other emissions 
Chemistry) CBM-IV* 


Tiedtke (1989) 

Louis (1979) 

Landgraf and Crutzen (1998) 
Krol and van Weele (1997) 
Wesely (1989) 

Guelle et al. (1998) 

Hao and Liu (1994) 

Price and Rind (1992) 

Levy et al. (1996) 

Olivier et al. (1996) 

Gery et al. (1989) 


Dentener et al. (1999) 
Dentener et al. (1999) 
Dentener et al. (1999) 
Lelieveld and Dentener (2000) 


Ganzeveld et al. (1998) 

Jeuken (2000) 

Marufu et al. (2000) 

Meijer et al. (2001) 

Lelieveld and Dentener (2000) 

Houweling et al. (1998) 


Houweling et al. (1998) 
a This column refers to previous work with the TM3 model. 
b This process is discussed in Section 2. 


archived diffusion and entrainment/detrainment coefficients 
from the ECMWF model using a modified scheme by Louis 
(1979) and Tiedtke (1989) scheme respectively. 

All the input fields are extrapolated from the ECMWF 
archived resolution (spectral, T106) to the TM3 model res- 
olution. The results presented in Section 5 and 6 were 
obtained using the multi-year (1979-1993) simulation by 
Lelieveld and Dentener (2000), which was also used in Pe- 
ters et al. (2001). The configuration of the TM3 model in the 
sensitivity studies (Section 8 ) is identical, and pertain to the 
year 1992. However, the model currently runs on a new SGI- 
3800 massive parallel system from the Dutch National Su- 
percomputing Facility (NCF). A short list of processes in the 
model and references to more elaborate discussion of these 
processes is given in Table 1. We will limit ourselves to a 
description of the processes explicitly relevant for this study, 
being the emissions from biomass burning, the photochem- 
istry of ozone and its precursors, and the lightning NO x pa- 
rameterization. 

2.1 Biomass burning emissions 

The emissions of biomass burning in the model were origi- 
nally constructed from two datasets. The EDGAR emission 
database (Olivier et al., 1996) prescribes the total amount of 
CO, HCHO, NO x , and NMHC emitted yearly from anthro- 
pogenic and natural sources. This database is based on es- 
timates for the year 1990 that were obtained from various 
sources on a resolution of 1° longitude x 1 ° latitude. For 
the temporal distribution of these emissions, the model re- 
lies on the Hao and Liu (1994) fire calendar, which contains 
monthly fractions of annual amounts of biomass burned on a 
5° longitude x 5° latitude resolution. For the sensitivity stud- 


ies presented in Section 8 , we have replaced the fire calendar 
of Hao and Liu (1994) with one based on fire counts from 
the AVHRR satellite (Flasse and Ceccato, 1996), which con- 
tains 10-day composites of global fire activity for 1992-1993. 
The observations were used to construct a time distribution 
function for each gridbox, by scaling the monthly total of 
counts by the yearly total of counts, yielding a fraction of 
fires occurring each month. Thus, the yearly totals remain 
the same, but the emissions maximize in the month with the 
largest number of fires. 

It is important to note that this approach is meant to pro- 
vide a simple alternative to the use of the Hao and Liu (1994) 
database, suitable for a model sensitivity analysis. In our al- 
gorithm, we do not account for the size of the fires that were 
observed; large fires and small fires contribute equally to our 
calendar, while in reality the amount of trace gases released 
to the atmosphere can differ strongly. Our algorithm will 
therefore overestimate the influence of these small fires and, 
since the total emissions are constant, underestimate the in- 
fluence of large fires. Also, our algorithm does not include 
the smoldering phase of the fires, which can continue for sev- 
eral weeks especially after forest fires [M.O. Andreae and P 
Merlet, Emission of trace gases and aerosols from biomass 
burning, Global Biogeochemical Cycles , accepted , 2000}. 
Due to the neglect of the smoldering phase, the timing of 
smoldering-related emissions (such as CO, NH 3 ) is biased 
towards the moment of detection of the fire, i.e. the flaming 
phase. Nevertheless, our alternative fire calendar suffices for 
a sensitivity analysis, whereas real satellite derived emission 
estimates require a more sophisticated approach (Pinty et al., 
2000) [B. Duncan, R. Martin , A. Staudt , R. Yevich and J. Lo- 
gan: Interannual and Seasonal Variability of Biomass Bum- 
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ing Emissions Constrained by Remote -Sense d Ob sensations, 
manuscript in preparation , 2001]. 

22 Photochemistry of ozone 

Photochemistry in the TM3 model is based on the CBM-4 
mechanism of Gery et al. (1989). It includes 38 trace gases 
and 110 gas phase reactions. The CBM-4 mechanism was 
initially intended to simulate ozone photochemistry over ur- 
ban areas with high NO x abundances. Therefore, Houwel- 
ing et al. (1998) modified the scheme such that low NO x 
conditions are better represented. This has improved ozone 
photochemistry over more remote areas, where NO x concen- 
trations are much lower. To reduce the number of reactions 
and species, the CBM-4 mechanism uses a form of struc- 
tural lumping, where trace gases of common origin with sim- 
ilar structures, lifetimes and reactive properties are treated as 
one species, which has the common characteristics of all the 
species it is meant to represent. This concept and the conse- 
quences for the photochemistry are described in more detail 
in the work of Gery et al. (1989), Duncan and Chameides 
(1998), Houweling et al. (1998) and Poschl et al. (2000). 

Precursors for ozone production that are emitted from 
biomass burning include NO x , CO, formaldehyde and ac- 
etaldehyde. Other NMHC compounds that are known to 
originate from biomass burning (e.g. acetone, acetonitrile, 
ethene) are represented through the emissions of several 
lumped species (parafinic and olefinic NMHCs, methylgly- 
oxal). The scheme includes the photochemistry of peroxy- 
acetyl nitrate (PAN), which is an important secondary prod- 
uct from the oxidation of acetaldehyde and contributes sig- 
nificantly to the free tropospheric nitrogen budget (Singh 
et al., 1995). Dry deposition is an important sink for NO x , 
HNO3, PAN, formaldehyde and acetaldehyde, while wet de- 
position through scavenging by clouds removes predomi- 
nantly HNO3, formaldehyde and H2O2. The chemical mech- 
anism includes the heterogeneous conversion of NO x to 
HNO3 on sulfate aerosols according to Dentener and Crutzen 
(1993), as well as a description of sulfur chemistry (Jeuken, 
2000 ). 

2.3 Lightning NO x emissions 

The lightning NO x emissions in the model are based on the 
widely used parameterization of Price and Rind (1992). The 
spatial distribution of lightning, as well as the intensity, is 
calculated using the heights of cloud tops in systems of deep 
convection. To discriminate between the efficiency of light- 
ning NO x production over the oceans relative to land, we 
have assumed a ratio of 1:10 in production efficiency, fol- 
lowing Levy et al. (1996). The lightning NO x production 
efficiency in cloud- to-ground discharges has been scaled to 
account for the higher efficiency in high latitudes (Macker- 
ras and Darveniza, 1994). The NO x emissions are uniformly 
distributed between cloud base and cloud top, proportional to 
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air-density. By scaling the lightning NO x emissions with me- 
teorological parameters, interannual variability in the model 
meteorology is also reflected in interannual variability in the 
lightning NO x emissions. The total annual simulated light- 
ning NO x emissions over the period 1979-1993 are 4.6 ± 0.3 
TgNyr- 1 . 

In our sensitivity experiments, we have used the verti- 
cal lightning NO x distribution profiles as given in Pickering 
etal. (1998). In this work, the vertical distribution depends 
on the location (tropics or extra- tropics), and surface charac- 
teristics (land or sea). Also, we have implemented a lightning 
NO x parameterization developed by Meijer etal. (2001). 
This follows the principles of Price and Rind (1992), how- 
ever, the occurrence and intensity of lightning NO x emis- 
sions is coupled to the convective precipitation in the model. 
This was shown to correlate well with space-based obser- 
vations of lightning flash rates over Europe, and was tested 
against NO x observations in the extra- tropics. The main dif- 
ference between the two approaches is the fraction of the 
emissions occurring in the tropics. In the new parameteriza- 
tion, this fraction decreases from 80% to 70%, which agrees 
better with independent estimates (Levy etal., 1996; Price 
and Rind, 1992). 

2.4 Labeled tracers 

To attribute photochemically produced ozone to the different 
processes in our model, we applied a tracer labeling tech- 
nique. The assumption that ozone formation in the tropo- 
sphere is NO x -controlled (Crutzen, 1988) implies that NO- 
to-NC>2 conversion by peroxy-radicals can be used to diag- 
nose ozone formation. Thereto, we have introduced labeled 
tracers of NO x , the most important NO y species, and ozone. 
The NO x tracers are labeled according to their source cate- 
gory, representing emissions from lightning, biomass burn- 
ing, nonagricultural soils, and industrial sources. When a la- 
beled NO molecule is converted to NO2 by a peroxy radical, 
the subsequent photolysis of NO 2 yields an ozone molecule 
labeled accordingly. Furthermore, tropospheric ozone intro- 
duced through stratosphere-troposphere exchange (STE) is 
labeled by tagging molecules once they pass the 100 hPa 
level. The labeled tracers participate in all other processes 
open to their regular, non-labeled counterparts. This ap- 
proach was tested and discussed in Lelieveld and Dentener 
( 2000 ). 

The sum of all the labeled tracers represents ~90% of the 
total ozone concentration in the tropical troposphere. The 
residual of ~ 10% that we cannot reproduce in the tracer la- 
beling is simply a result of the chosen approach. The la- 
beled ozone tracers in the model are a diagnostic tool, and 
not an exact solution of the differential equations govern- 
ing ozone photochemistry. Although many reaction path- 
ways, including cycling of NO x through several reservoir 
species, are labeled, not all possible routes and recombina- 
tions are taken into account. For instance, cycling of NO2 
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and NO3 through N2O5 can yield back the original products, 
or HNO3. The many combinations in which this is possible 
given five marked tracers (plus one unmarked) is large, and 
had to be simplified. Also, the production of marked ozone 
from NO+RO2 reactions is calculated simply from the reac- 
tion fluxes, which can also lead to errors as was discussed 
extensively in the original work of Johnston and Kinnison 
(1998), and comments on that work Poschl et al. (2000). For 
the purposes in this paper the labeled tracer approach, al- 
though clearly not exact, gives very useful information and 
the ability to diagnose the ozone budget qualitatively. 

In the calculations presented in the next sections, the local 
tropopause height is calculated in the same way as in Pe- 
ters etal. (2001). We define the tropical tropopause as the 
level where the temperature lapse rate changes sign. This 
gave better agreement with the observed tropopause heights 
(see Peters etal. (2001)) than gradient based methods such 
as Craig (1965). Due to the coarse vertical resolution, our 
model underestimates tropopause heights in the tropics by 
about 20 hPa, which translates to a 1 .6 DU underestimate of 
tropospheric ozone columns assuming a mixing ratio of 100 
hPa at the tropopause. 

3 Observations 

We use the NIMBUS-7 TTOC record (1979-1992) to detect 
systematic flaws in our model and to gain more insight in the 
sensitivity of our model calculated ozone. The observations 
were derived from the TOMS total ozone observations with 
the Modified-Residual (MR) method (Hudson and Thomp- 
son, 1998). 

The TTOC record provides twice monthly fields of TTOC, 
covering all the tropical latitudes where a zonal wave one pat- 
tern in total column ozone is observed. The horizontal res- 
olution of the dataset is 2° longitude xl° latitude (230x115 
km at the equator), and through averaging over several or- 
bital passes of the satellite, a picture of the full zonal distri- 
bution of TTOC is obtained in several days. Through aver- 
aging the observations over 13-15 day periods, errors asso- 
ciated with the scan angle, as well as instrument noise are 
filtered out. Thus, the product is not suited to study day-to- 
day variability. The data is available on the internet through 
http://metosrv2.umd.edu/~tropo/. 

Uncertainties in the observations of TTOC are related to 
the assumed tropopause height as well as the insensitivity 
of the TOMS instrument to tropospheric ozone near the sur- 
face. The instrument is progressively insensitive to the signal 
of ozone originating from altitudes below the stratospheric 
ozone maximum. The algorithms that were developed to cir- 
cumvent these problems are explained in detail in Hudson 
and Thompson (1998). In the MR retrieval algorithm, sonde 
observations from selected stations in the SH are used to ob- 
tain a-priori information on the zonal distribution of tropo- 
spheric ozone, which is thereafter projected onto each grid 


point of the total column ozone observations to yield the 
TTOC. The averaging of several sondes to obtain a half- 
month average leads to a standard deviation of ±5.3 DU. 
Therefore, the a-priori zonal ozone distribution cannot be 
fixed more accurately than this sonde standard deviation, and 
5.3 DU is thus the lower limit for the precision in TTOC. 
This value is comparable to the TOMS instrument accuracy. 
An independent comparison of TTOC with four ozone sonde 
stations in the tropics revealed a RMS deviation of 5.1 DU 
between TTOC and sondes (Thompson and Hudson, 1999). 
This value is expected to increase with distance from the 
sonde stations used in the algorithm, and is probably larger 
in the NH than in the SH, since the underlying zonal wave 
one pattern used in the retrieval was based on SH ozone son- 
des. However, since there is no objective way to define the 
deviations at other locations, the overall standard deviation 
of the TTOC product is given by cr=5.3 DU over the whole 
domain. 

The Nimbus-7 TTOC record (1979-1992) has been evalu- 
ated in two ways with ozonesonde data. Note that these com- 
parisons are mostly performed in the Southern Hemisphere 
tropics [A. M. Thompson et al , The 1998-2000 SHADOZ 
( Southern Hemisphere ADditional OZonesondes) tropical 
ozone climatology. 2. Stratospheric and tropospheric ozone 
variability and the zonal wave-one , J. Geophys. Res., sub- 
mitted, 2002, from now on: AM. Thompson et al, submit- 
ted manuscript, (2002), since the available Northern Hemi- 
sphere tropical data are restricted to campaign data (Thomp- 
son etal., 2000) and to stations at Suriname (5N, 55W) and 
Malaysia [S. Yonemura et al.. Tropospheric ozone climatol- 
ogy over peninsular Malaysia from 1992 to 1999, J. Geo- 
phys. Res., in press , 2002). First, the limited ozonesonde 
data for the 1979-1992 period agreed with the TTOC satellite 
product within 6-7 DU on average (Thompson and Hudson, 
1999). Second, a statistically averaged annual dataset was 
prepared from Nimbus-7 TTOC and compared to the 1998- 
2000 SHADOZ ozonesonde record in four regions: tropical 
Atlantic, central Pacific, eastern Pacific and Indonesian mar- 
itime continent. Excellent agreement in TTOC magnitude 
and seasonality was found in the central Pacific TTOC com- 
pared to sondes at Samoa (14S, 17 1W). The TTOC season- 
ality over the eastern Pacific and Indonesian maritime conti- 
nent compared well to ozonesonde seasonal patterns at San 
Cristobal (IS, 89.6W) and Watukosek, Java (7.5S, 1 13E), re- 
spectively, but tended to be too high in magnitude. Over 
the Atlantic, comparison of the TTOC annual mean with 
ozonesondes over Ascension Island (6S, 15W) suggests an 
overestimate of 5-8 DU on average during the period from 
July-October 1998-2000. 

It should be noted that differences among several retrieval 
algorithms of TTOC (e.g. Ziemke etal. (1999), Fishman and 
Larsen (1987)) can be quite large, also compared to the re- 
ported standard deviation of 5.3 DU. No comprehensive in- 
tercomparison has been published so far, although Martin et 
al., [R. V. Martin et al. Global model analysis of TOMS and 
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Fig. 1 . Observed (top) and mode] calculated (bottom) distribution of TTOC in September 1992. Note the widespread underestimate of TTOC 
over the tropical Atlantic Ocean and adjacent continents. 


in -situ obse nations of tropical tropospheric ozone , J. Geo - 
phys. Res., submitted, 2002 ] have noted a tendency for trop- 
ical Atlantic ozone from TOMS-based data to run high rela- 
tive to in-situ measurements, particularly in the second half 
of the year. Our experiences with working with the Ziemke 
et al. (1999) data suggest that the MR method generally gives 
higher values of TTOC, and several regions are suspected to 
experience unwanted aerosol influences in the MR method 
(A. Frolov, personal communication ). A detailed quantifica- 
tion of the uncertainties of individual locations in the MR 
method is beyond the scope of this study, as is an intercom- 
parison between the different TTOC datasets. More insight 
in the value and represen tivity of these datasets can be ob- 
tained through actively using them, as we do in this study. 
For now, we therefore take the MR-TTOC data ’as is’, and 
are cautious not to overinterpret our results. 

Alternatively, ozone sondes and in situ measurements 
from campaigns such as TRACE- A, PEM, NOXAR, or 
MOZAIC could be used to test our models performance in 
the tropics. The data acquired in these campaigns generally 
include detailed vertical information on several trace species 
and meteorological parameters. However, there is an impor- 
tant trade off to consider that can be seen as a choice between 
the amount of detail of information and the areal and tempo- 
ral coverage. The individual campaigns usually cover peri- 
ods of several weeks in a specific year, and focus on a limited 
geographical region. Especially near the surface, they are of- 
ten influenced by local pollution or particular small scale me- 
teorological conditions (e.g. sea-breeze, orographic flows). 
This is especially disadvantageous when comparing to global 
models that generally apply coarse resolutions and parame- 


terizations not designed to capture these small scales. Ozone 
sondes and in situ measurements are much better suited to 
study smaller temporal and spatial scales, or the effect of a 
single process (e.g. biomass burning or convection). These 
data would therefore form an ideal basis to improve possible 
model shortcomings, which is the next step in our research. 
However, in this work we take advantage of the global, multi- 
year TTOC record to firstly isolate the regions, seasons, and 
processes that require improvement. 

4 Measure of error 

Figure 1 (top) shows the observed TTOC distribution in 
September 1992, when biomass burning emissions were par- 
ticularly strong. The TOMS-based TTOC is in excellent 
agreement with in-situ observations (Thompson etal., 1996; 
Thompson and Hudson, 1999). Ozonesondes at Ascen- 
sion Island and the near-coast locations of Brazzaville (4S, 
15W) and Natal (6S, 35W), averaged 50 DU in September 
1992. On aircraft flights out of Ascension during October 
1992, TTOC from lidar measurements ranged from 45-55 
DU (Browell et al., 1996). 

The South Atlantic ozone maximum clearly stands out 
with values >50DU, gradually decreasing with distance from 
the Atlantic, reaching minimum values between 10 and 20 
DU over the Pacific Ocean. The comparison with model sim- 
ulation of TTOC for the same year and month in Fig. 1 (bot- 
tom) shows distinct differences. The absolute values of the 
TTOC often differ by more than the one sigma (5.3 DU) stan- 
dard deviation of the measurements, and many smaller scale 
features are not reproduced. The discrepancies are largest 
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over the Atlantic Ocean, where a systematic mismatch is seen 
between observations and model calculations. In contrast, 
the model reproduces the TOMS observations quite well over 
the Pacific Ocean. This is consistent with earlier results (Pe- 
ters etal., 2001), based on a comparison with ozone sondes 
from Samoa, Ascension and Natal. 

Beside this ’Atlantic mismatch’, there are two other im- 
portant features that appear from Fig. 1. Firstly, the amount 
of detail in the observations is larger. This can be partly at- 
tributed to the higher resolution of the observations, but it 
is also caused by a lack of variation in the model’s monthly 
emission data. To account for the different resolutions of 
the datasets, we interpolated the model results to the 2°xl° 
degree resolution of the TTOC observations and we applied 
a zonal running average to the TTOC observations, with a 
window size equal to the zonal model grid-size (5° degrees). 
This filters out strong local features in the observations at 
resolutions smaller than the model grid, which cannot be re- 
produced by the model. 

Secondly, it is important to note that variations in the 
meridional direction are small compared to the large zonal 
variation. Since the area depicted in Fig. 1 encompasses 
only 5-6 model grid boxes in the meridional direction (appr. 
15° degrees), the model can not reproduce fine scale struc- 
tures in the meridional direction. Moreover, the zonal TTOC 
distribution is typically a feature on large (planetary) scale, 
while meriodional variations within the tropics are on a much 
smaller, almost local scale. Therefore, the next step in our 
comparison is to take the meridional average of the respec- 
tive datasets at each longitude, which yields the zonal distri- 
bution of TTOC. 

4.1 Zonal TTOC distribution 

Figure 2 shows the derived zonal TTOC distribution in 
September 1992 for both datasets. It has already been 
pointed out that the MR-TTOC for September 1992 agrees 
well with the available sonde data near the Atlantic (location 
of peak in zonal maximum. Fig. 2). More comprehensive 
ozonesonde data, collected in 1998-2000 from 10 stations 
across the entire tropical zone, confirm the wave-one pattern 
(A. M. Thompson et ai f The 1998-2000 SHADOZ (Southern 
Hemisphere ADditional OZone sondes) tropical ozone Cli- 
matology ; L Comparison with TOMS and ground-based 
measurements, J. Geophys. Res., in press, 2002, from now 
on: AM. Thompson et al, in press, (2002)). A longitudinal 
cross-section of those observations, shows the wave to be as- 
sociated with elevated ozone (> 45 ppbv) concentrated in the 
free troposphere between 50W and 50E (AM. Thompson et 
al, submitted manuscript , 2002). The distribution and inten- 
sity of the higher-ozone feature, hence the magnitude of the 
wave, varies with time of year, but it is always present. The 
regions of enhanced ozone correspond to locations where 
ozonesondes and other observations implicate biomass burn- 
ing, lightning (Thompson etal., 1996), large-scale subsi- 
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Fig. 2. The meridionally averaged zonal TTOC distribution in 
model calculate4 (blue) and satellite retrieved (red) TTOC for 
September 1992. This distribution is used to define a measure of 
error according to Eq. (1). The shaded lines denote one-sigma stan- 
dard deviation (5.3 DU). Atlantic ozonesonde and aircraft ozone 
data (taken within longitudes 35W to 15E) during September and 
October 1992 fall within the shaded region. 


dence (Krishnamurti etal., 1996) and stratospheric intrusions 
(Baray etal., 2000) as possible ozone sources. 

The figure clearly shows the zonal extent of the Atlantic 
mismatch, characterized by model calculated TTOC values 
that are consistently lower than the observations over the 
90W to 90E region. Although this mismatch maximizes each 
year in September, other months show a similar problem. 
The TTOC distribution depicted in Fig. 2 is chosen as the 
focal point of our comparison. Conveniently, the TTOC re- 
trieval algorithm uses the meridional extent of the zonal wave 
one pattern to distinguish tropical air masses from extra- 
tropical air masses. Thus, this pattern is always present in 
the satellite observations, and points that do not qualify as 
’tropical’ in the retrieval algorithm are not included in the 
comparison to our model calculations. 

To make an objective comparison between the observed 
and modeled zonal TTOC distributions, a mathematical mea- 
sure of error is introduced, which consists of two compo- 
nents. First, we define c lf which is simply the correlation (r) 
between the observed and modeled zonal wave one pattern 
in Fig. 2. The parameter ci indicates how well the horizon- 
tal gradients in the zonal TTOC distribution are reproduced, 
and how well transitions in TTOC induced by the land- sea 
distribution for instance are captured. The shape of this pat- 
tern is influenced amongst other factors by the location of 
the Inter Tropical Convergence Zone (ITCZ), the presence 
of biomass burning plumes, or the occurrence of large scale 
convection and associated lightning. Increasing values of ei 
denote a better representation of these zonal gradients. It is 
important to note that ei is not influenced by the absolute 
value of TTOC, and its value will increase only when the 
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zonal gradients in the model calculations better fit the obser- 
vations. This is advantageous because it excludes the errors 
in the measurements and model that are introduced through 
uncertainties in the tropopause height. On the other hand, the 
parameter £1 cannot be used to assess the absolute fit between 
the observations and the calculations. 

Therefore, we also introduce a second parameter, €2, that 
is defined as the absolute deviation between observations and 
measurements, scaled by the standard deviation of the obser- 
vations: 


e 2 (<) = 


1 

N lon 


£ 


TOMS{t) - TM3(t) 

&TOMS 


( 1 ) 


Where TOMS(t) is the obs erved zo nal distribution of 
TTOC for each month (t), and TM3(t) is the model calcu- 
lated zonal distribution of TTOC for each month. The value 
of €2 will decrease as the discrepancy between observations 
and calculations decreases. 

The number of longitudes used in the summation can vary 
depending on the scale we want to look at. For instance, 
in the analysis of the variability (Section 5), we analyze 14- 
year time series of 62 at each location (N/on=l, t=168), while 
for the analysis of our sensitivity runs (Section 8) we sub- 
divide the zonal domain into only two regions; a Pacific 
region (N ion =90 , t=168) and an Atlantic region (Ni on =90, 
t=168). Together, e\ and €2 enable us to quickly and objec- 
tively quantify the agreement with the TTOC observations, 
and intercompare the results of different model simulations. 


4.2 Algorithm test 

In order to assess the robustness of our algorithm, and to get 
an indication of the sensitivity of e\ and €2, we have applied 
our measure of error to a series of test distributions of ozone. 
Thereto, an ensemble of one hundred 3-dimensional ozone 
distributions was created, in which each gridpoint of an ac- 
tual model calculated ozone distribution was randomly per- 
turbed by an average of ± 10 %. These random perturbations 
represent the effects of for instance anomalous winds, con- 
vective venting, or unexpected emissions of ozone precursors 
that could disturb the ozone concentration within a gridbox 
in a non-systematic way, that is nevertheless substantial in 
magnitude. Next, we have integrated the columns, averaged 
to make the zonal pattern, and then calculated the yearly val- 
ues of ei and €2 for each of the one hundred perturbed dis- 
tributions. The distribution of and 62 in this ensemble re- 
veals the spread that can be expected from relatively large, 
but random perturbations. Any change in £1 and £2 that is 
significantly larger (more than lcr) than this, should then be 
seen as a systematical change of our model calculated ozone 
field. 

From our one hundred member ensemble, we have derived 
a one-sigma standard deviation of ± 0.00 1 for 61, and ±0.01 
for £2. We can thus say that the small differences we report 
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Fig* 3. Hovmoller plot of €2 as a function of longitude for the pe- 
riod 1979-1992. A value of 1 denotes a difference of 5.3 DU (1-cr) 
between model and observations. The +/- signs were added to the 
figure to signify model underestimates or overestimates repectively. 


in Section 6 and Table 2 represent a systematic difference 
between the calculated ozone fields from different years, and 
under different scenarios. 


5 Comparison of TTOC 

As a first step, we calculate the values of £2 for each longi- 
tude and each month (t=168) of our multi-year simulation. 
The result is shown in Fig. 3, where a negative sign indicates 
that the model calculated values are higher than the observa- 
tions. It shows that the values of £2 range from 0-2 a (appr. 
0-10 DU) over the tropical Pacific, to more than 3 a (appr. 
16 DU) over the tropical Atlantic. The largest underesti- 
mates occur during August through November. The model 
calculated values of TTOC are generally lower than the ob- 
servations. The only time the model overpredicts the TTOC 
values is in Dec-Jan-Feb over the Pacific Ocean. 

It is clear from Fig. 3 that the Atlantic mismatch is the 
most serious problem in our comparison, and that it is an an- 
nually recurring feature. In contrast, high values of £2 over 
the Pacific Ocean appear for several years, but not in all. To 
separate these systematic differences from the ones that oc- 
cur occasionally, we have decomposed each time series of 
£2 into a seasonal cycle, and a component of interannual 
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Year of TOMS observations 


Fig. 5. Correlation coefficient (e , ) between the zonal TTOC distribution in observations and simulations. For each year of TTOC observations 
(1979-1992) correlations were calculated with each year of a multi-year simulation with our CTM (1979-1992). The bars denote the range 
of correlations, the red diamonds are the averages. This value is representative for a simulation with ’climatological’ meteorology. The green 
diamonds signify the correlation when the correct meteorology of the year of observations is applied to the model. 


variability. For a discussion of the interannual variability in 
TTOC we refer to Peters etal. (2001). The seasonal cycle 
of €2 is shown in Fig. 4. The figure confirms that most of 
the Atlantic mismatch is caused by seasonal differences. On 
average, the seasonal cycle contains 75% of the variance of 
the €2 time series over this region. In contrast, the seasonal 
cycle of 62 over the Pacific Ocean explains less than 50% 
of the original variance. If we consider the relative errors, 
both deviations translate to ca 30% of the observed TTOC. 
From this, we can conclude that the most important model 
improvements should focus on the seasonality of TTOC over 
the Atlantic Ocean. Possibly, the Pacific region would profit 
from such improvements as well, since the two regions are 
coupled through transport. 


6 Influence of transport 

The coupling between the regions through transport is not 
the same in each year due to the occurrence of anomalous 
meteorological conditions such as those induced by ENSO. 
Through the use of ECMWF wind fields for each year, we 
have included this transport variability, and in Peters etal. 
(2001) showed that the effect of ENSO on TTOC is present 
in our multi-year simulation. In this section, we will apply 
our measure of error to quantify the influence of transport on 
our results. 


We compare each year of TTOC observations with each 
year of our multi-year simulation. In our multi-year sim- 
ulation, we have 15 years of calculated TTOC, which dif- 
fer in meteorology (including lightning), stratospheric upper 
boundary conditions for ozone, and the anthropogenic emis- 
sions which include trends over the 1979-1993 period. In 
this section, we use yearly averaged values for 1979-1992 
of the monthly calculated 6 i to assess the performance of 
our model. Since ei is independent of the absolute values 
of TTOC, the influence of increasing emissions in the multi- 
year simulation is largely excluded in this comparison. Thus, 
any year-to-year difference in ei is due to changes in meteo- 
rology (including changing lightning patterns). 

For each year of the TOMS observations (1979-1992, 
n^M), we have calculated 14 values of ci; one value for each 
year of simulations (1979-1992, m=14). These individual 
values of €i denote the success in simulating the observed 
zonal distribution of TTOC, given the specific meteorolog- 
ical conditions of each of the 14 simulated years. The av- 
erage 6 i of these individual years would represent the value 
when a random year of meteorology is applied, while there is 
one value of ei that corresponds to the correct meteorology 
matching the year of TTOC observations. 

Figure 5 shows the multi-year average value (open dia- 
monds) and range of interannual variation (bars) for each of 
the 14 years TTOC observations, as well as the values of 
ei for the model year corresponding to the year of observa- 
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Fig. 4, 14- Year average seasonal cycle of €2 as a function of lon- 
gitude. A value of 1 denotes a difference of 5.3 DU (1-cr) between 
model and observations. The +/- signs were added to the figure to 
signify model underestimates or overestimates repectively. 


tions (filled diamonds). The simulations performed with ac- 
tual meteorology consistently perform better than average. In 
fact, the simulation with actual meteorology generally ranks 
higher than any other year. This is due to the fact that me- 
teorological conditions determine to a large extent the loca- 
tion of local ozone maxima, as well as the zonal gradients 
in TTOC. Although the differences between a correlation of 
0.7 and 0.84 might seem small, the fact that these differences 
show up in a yearly and meridionally averaged TTOC dis- 
tribution indicates that the individual differences are in fact 
quite large, and certainly significant given the sensitivity of 
e\ (Section 4.2). 

From the figure one can also see that the 14-year aver- 
age value of ei varies from year-to-year, as does the range. 
This shows that some years of observations are more diffi- 
cult to model than others presumably because their specific 
conditions are more difficult to capture. This is obvious for 
1982/83, when El Nino conditions prevailed. While the sim- 
ulation with actual meteorology performs well, some other 
years show a poor correlation. This indicates that in such a 
year, the zonal TTOC distribution can not be modeled very 
well if one does not consider the specific meteorological con- 
ditions during ENSO. The fact that the maximum of €1 for 
1983 is the lowest in the whole series shows that other pro- 
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Fig. 6. Absolute contribution of several processes to the modeled 
TTOC distribution. Processes include photochemical production 
of ozone following NO x emissions from soils, lightning, biomass 
burning and industrial sources, as well as ozone from stratosphere- 
troposphere exchange. Values for September 1992. 


cesses than transport that change during ENSO also influ- 
ence the shape of the zonal TTOC distribution. Such pro- 
cesses include the occurrence of wildfires following anoma- 
lous droughts. Although this was shown to have a large im- 
pact on TTOC (Chandra et al., 1998; Thompson et al., 2001), 
it was not included in our multi-year model simulations. 
Adding this variability in emissions to the model is a dif- 
ficult challenge, for which tools are currently being devel- 
oped [ B . Duncan ; R. Martin , A. Staudt, R. Yevich and J. Lo- 
gan: Interannual and Seasonal Variability of Biomass Burn- 
ing Emissions Constrained by Remote -Sensed Observations, 
manuscript in preparation, 2001). 


1 Role of photochemistry 

Beside the previously mentioned wildfires, other sources 
contribute to the photochemical production of ozone in the 
troposphere, and thus partly determine the zonal distribu- 
tion of TTOC as depicted in Fig. 2. By using the labeled 
tracers (see Section 2.4), we have calculated the contribu- 
tion of several tropospheric ozone sources to the modeled 
TTOC distribution. The processes that were thus separated 
are photochemical ozone production from lightning, indus- 
trial, biomass burning, and soil NO x emissions, as well as 
the contribution from the stratospheric influx of ozone. Fig- 
ure 6 shows the contribution of these processes to the mod- 
eled zonal TTOC distribution in September 1992. 
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Ozone production from industrial NO x emissions con- 
tributes several DU to the modeled ozone concentration. 
There is only very little difference in zonal direction, be- 
cause the most important source is located in the NH extra- 
tropics. Transport in longitudinal direction by the power- 
ful jet streams of the NH and SH occurs more rapidly than 
meridional transport to the tropics and across the ITCZ, lead- 
ing to a largely zonally invariant distribution in the tropics. 
Lightning NO x emissions in the tropics are stronger over 
the continents than over the oceans. Yet, the distribution of 
TTOC from lightning NO x does not show much structure 
in zonal direction, because the NO x enriched air is brought 
into the middle and upper free troposphere, where trans- 
port is fast, and the lifetime of ozone relatively long. The 
monthly mean signal of ozone from lightning NO x is there- 
fore seen throughout the tropics, adding 5-8 DLT of ozone 
to most columns. The strongest influence is seen over the 
Atlantic Ocean and South American content, as is expected 
during this time of year. 

An almost equal contribution to TTOC comes from 
stratosphere-troposphere exchange. However, most of this 
ozone is located at the edges of the domain in our analy- 
sis, where the subtropical jet sometimes penetrates far equa- 
torwards. This happens at locations where the jetstream is 
disturbed most frequently; near surface orography and land- 
sea barriers. Thus, most of the stratospheric ozone encoun- 
tered in the tropical troposphere in our model has crossed 
the tropopause in the sub-tropics, and is subsequently trans- 
ported equatorwards from these locations where disturbances 
in the jetstream occur. 

The strongest contributors to the ’wave one’ character of 
tropical TTOC are soil NO x emissions and biomass burning 
NO x emissions. Both these sources are terrestrial surface 
sources (50% of soil emissions in the tropics), which maxi- 
mize over South America and Africa. The ozone produced 
from these NO x emissions is transported more slowly in the 
BL, lower, and middle free troposphere, and has a shorter 
lifetime there. It therefore slowly accumulates over the At- 
lantic, but exerts only little influence on the more distant trop- 
ical Pacific region. 

The seasonality of ozone from biomass burning (not shown) 
is very similar to the seasonality of 62 in Fig. 4. In contrast, 
ozone produced from lightning NO x emissions displays its 
maximum in April-May-June, while ozone produced from 
soil emissions is almost constant throughout the year. This 
again suggests that biomass burning could be partly respon- 
sible for the observed Atlantic mismatch. To see which part 
of our treatment of biomass burning could be the cause of the 
seasonal ozone underestimate, we have performed a sensitiv- 
ity analysis. 
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Table Z Experiments conducted and their values of the parameter 
Ci or 62 * Experiments that meet the criteria described in Section 8 . 1 
are displayed in boldface, and discussed in the text. Improvements 
in €1 are achieved from a more than 0.001 increase in value, and 
improvements in €2 from a 0.01 decrease. 


Experiment 

Cl 

£2 .Pac 

62 ,Atl 

Base 

0.842 

0.77 

1.80 

Biomass 

burning 



2x NOx 

0.848 

0.71 

1.58 

2x CO 

0.843 

0.73 

1.71 

2x NMHC 

0.844 

0.76 

1.77 

2 x J(HN0 3 ) 

0.702 

0.82 

1.79 

Injection height 

0.853 

0.72 

1.52 

Alternative calendar 

0.852 

0.76 

1.74 

Lightning 



2x NO x lightning 

0.848 

0.52 

1.08 

5 x NOx oceans 

0.723 

0.62 

1.26 

Pickering profiles 

0.797 

0.69 

1.67 

Coupling to conv. prec. 

0.814 

0.82 

1.81 

land/ocean ratio 2: 1 

0.762 

0.70 

1.85 

Other 



5 DU offset 

0.842 

0.57 

1.01 


8 Sensitivity analysis 

Beside our treatment of biomass burning emissions, the sen- 
sitivity analysis includes simulations to investigate uncertain- 
ties in the lightning NO x emission parameterization. These 
two processes occur predominantly in the tropics, in contrast 
to industrial emissions, soil emissions and STE. Although 
the seasonality and zonal distribution of lightning NO x emis- 
sions do not match the pattern of £2 (Fig 4), its large contribu- 
tion to the ozone budget and large uncertainties in magnitude 
and distribution justifies a more detailed look at its parame- 
terization. 

8.1 Strategy 

To ensure an objective intercomparison between different 
model simulations, we now define £2 for two different 
regions. The first region encompasses all longitudes between 
90W and 90E (^^=90, ’Atlantic’), and the other region 
encompasses the longitudes outside this area (180W-90W 
and 90E-180E, N ion =90, ’Pacific 1 ). Since the good perfor- 
mance of the model over the Pacific Ocean was determined 
independently by comparison with ozone sondes (Peters 
etal., 2001), our strategy to determine whether a simulation 
constitutes an improvement over our base run is: 1) ensure 
that €2 for the Pacific fit does not increase; 2) achieve a 
decrease of €2 for the Atlantic mismatch, and 3) try to 
increase the value of €\. Only when these three criteria are 
met, a systematic improvement of modeled TTOC has been 
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achieved. This has turned out to be a very strict constraint 
due to the large number of observations and the scale of the 
phenomenon. All the sensitivity simulations pertain to the 
year 1992. 

Table 2 presents an overview of all the sensitivity simu- 
lations that were performed, and their yearly average values 
of €i and 62 The base simulation is the model as used in Pe- 
ters etal. ( 2001 ) as well as in the previous sections of this 
work. The simulations that meet the criteria described above 
are highlighted in boldface. These simulations are discussed 
briefly in this section. We repeat that the aim of these simula- 
tions was to prioritize directions for further research. There- 
fore, the analysis and discussion of these simulations will be 
only brief, and simulations that did not qualify as improve- 
ments will not be addressed. 

Three of the possible improvements are related to biomass 
burning. Simulations were performed with an increased in- 
jection height of biomass burning products, double NO x 
emissions from biomass burning, and an alternative fire cal- 
endar, respectively. These experiments increase the amount 
of NO x in the free troposphere over the tropical Atlantic, and 
thereby increase the photochemical production of ozone. The 
amount of NO x reaching the free troposphere is determined 
by the complex interaction between emissions, deposition to 
the surface, oxidation by OH, and convective venting from 
the planetary boundary layer (PBL). 

The efficiency by which an NO x molecule contributes to 
ozone formation depends to a large extent on its lifetime, i.e., 
its ability to reach the free troposphere and avoid dry de- 
position or oxidation to HNO 3 . Due to the dry conditions 
prevailing during biomass burning events, strong convective 
systems can be assumed to be absent near large fires. More- 
over, NO x has a lifetime of only several days in the PBL, 
and in its oxidized form (HN0 3 ) surface deposition is even 
faster and rainout during convective transport occurs quickly. 
Therefore, the interaction between meteorology and emis- 
sions largely determines whether a NO x molecule will reach 
the free troposphere, or will be oxidized and/or deposited in 
the PBL. Our three successful experiments influence this in- 
teraction in three different ways. 

Firstly, we reduce rapid deposition of NO x in the PBL 
and increase the concentration in the free troposphere by 
distributing the biomass burning NO x emissions over the 
lower 6 km of the troposphere. Obviously this is arbitrary, 
but it shows that an improvement would be realized if more 
biomass burning NO x were transported to the free tropo- 
sphere in our model. This could indicate a lack of trans- 
port of biomass burning products from the planetary bound- 
ary layer (PBL), which is plausible since this parameteriza- 
tion depends on 6 -bourly updated diffusion parameters, that 
can not adequately capture the diurnal growth and decline of 
the PBL. Using Rn 222 , Dentener et al. (1999) already noticed 
that our model appears to underestimate convective venting 
of the PBL over the African continent. It could also indi- 
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cate that biomass burning smoke plumes may reach the free 
troposphere on their own buoyancy, which is not included in 
our model. Finally, Yienger and Levy (1995) suggested that 
post-burning soil emissions of NO x could contribute to the 
emission budget. These emissions could easily reach the free 
troposphere under unstable or convective conditions occur- 
ing after the fire activity has ceased. 

In addition, the interaction between transport and biomass 
burning in the model could he improved through implemen- 
tation of a more realistic fire calendar. As a result, the loca- 
tion and occurrence of convection would become more con- 
sistent with the emissions. The convective vertical exchange 
of NO x decreases in July and August, and increases by 15% 
in September and October. This leads to more ozone produc- 
tion in the free troposphere during these months. This agrees 
well with the findings of Cooke etal. (1996) and Galanter 
et al. ( 2000 ), who applied more detailed alternative fire cal- 
endars to point out shortcomings in the inventory of Hao and 
Liu (1994). The fact that ei increases significantly indicates 
that the interaction between biomass burning and convection 
has a strong influence on the zonal TTOC gradients. 

The third important parameter in the description of 
biomass burning is the emission magnitude. A doubling of 
the emissions of NO x from biomass burning also increases 
the flux of NO x from the planetary boundary layer. Table 3 
shows the consequences for the tropospheric ozone budget. 
The increase in production of ozone is largely balanced by an 
increase in destruction, as well as by an increase in deposi- 
tion and an increase in advection of ozone out of the Atlantic 
region. The net effect of a doubling of the NO x emissions 
is therefore a 2 Tg increase in ozone, which is only 5% of 
the tropospheric burden. The effect on the TTOC columns 
is similar (appr. 2 DU). The effect of increased emissions 
is more pronounced when it is coupled to the alternative fire 
calendar. In fact, this simulation shows the highest sensitiv- 
ity. This suggests that the correct representation of the occur- 
rence and magnitude of biomass burning is an important pre- 
requisite to model TTOC, in agreement with Galanter et al. 
( 2000 ). 

An alternative method to increase the concentration of 
NO x in the free troposphere is through lightning NO x emis- 
sions. Our current emission magnitude of 4.6 TgN yr 1 is 
in the lower part of the range of Inter Governmental Panel on 
Climate Change (IP CC) estimates (4-20 TgN yr“ J (Prather 
etal., 2001 )), and a doubling would still be within the un- 
certainty limits (although at die high end of currently used 
values in global models). Table 2 shows that this would yield 
an improvement in our comparison to TTOC. Clearly, the im- 
pact of this parameter is very large, and an increase in trop- 
ical lightning NO x emissions would give better agreement 
with TTOC observations. However, this would also strongly 
affect the concentrations of ozone in the extra- tropics where 
—20% of the emissions occur in our model. This would lead 
to a large overestimate of ozone over Europe, the United 
States, and Asia, where ozone concentrations were previ- 
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Table 3. Budgets of ozone in the base run, and in a simulation 
with double NO x emissions from biomass burning. The last col- 


umn shows the difference between both simulations . Numbers are 
tropospheric integrals over the tropical Atlantic region and adjacent 
continents [80W-40E,12S-12N] n for the year 1992. 

Base (A) double NO x (B) B-A 

Chem. production 

1095 

1200 

+ 105 

NO+H0 2 (%) 

62 

62 

0 

NO+CH 3 O 2 (%) 

16 

16 

0 

NO+XO 2 w 

17 

17 

0 

NO+CH 3 C 002 (%) 

5 

5 

0 

Chem. destruction 

951 

1019 

+68 

O 3 + hv (%) 

59 

58 

-1 

0 3 + OH (%) 

10 

10 

0 

0 3 + HO 2 (%) 

25 

26 

+1 

O 3 + HC (%)‘ 

6 

6 

0 

Net chemistry 

+ 144 

+181 

+37 

Advection 

26 

7 

-19 

Dry deposition 

166 

182 

+16 

STE C 

0 

0 

0 

Total 

44 

+6 

+2 

Tropospheric burden 

40 

42 

+2 


a Units are Tg0 3 yr“\ unless percentage explicitly stated. 
b HC denotes isoprene, ethylene and olefinic hydrocarbons. 
c direct STE in this region is negligible. 


ously modeled well (Lelieveld and Dentener, 2000). Al- 
though this good agreement could be a fortuitous result of 
several uncertain processes, we cannot neglect this impor- 
tant side effect of an increase in lightning NO x emissions. 
In the alternative parameterization of Meijer etal. (200 l) t 
the percentage of extra-tropical emissions would be even 
higher (30%), and the sensitivity over the Northern Hemi- 
sphere even larger. Alternatively, we adjusted the lightning 
NO x emissions through changing the lancVsea emission ra- 
tio, the emission magnitude over the oceans, and through in- 
clusion of vertical distribution profiles of emissions accord- 
ing to Pickering et al. (1998). However, none of these simula- 
tions met the criteria of a systematic improvement. For now, 
we will therefore only state that TTOC in the tropics is very 
sensitive to lightning NO x emissions, and further studies are 
urgently called for. 

9 Vertical ozone profiles 

To illustrate the results of our sensitivity experiments we 
present a comparison with selected vertical ozone profiles 
from several tropical stations. These stations are part of 
the SHADOZ program, which encompasses 1 1 ozone sonde 
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stations throughout the tropics. The details of the differ- 
ent stations and a comparison to TOMS data and ground- 
based measurements can be found in AM. Thompson et al., 
in press , 2002. 

Figure 7 shows six observed (black solid lines) vertical 
ozone profiles from selected stations and dates from the year 
1998, as well as several model calculated profiles for the 
same day in 1998. These modeled profiles show results from 
a standard model simulation (blue dashed lines), a simulation 
with the meteorology of 1997 applied (red dashed lines), and 
a simulation with doubled biomass burning NO x emissions 
and an alternative fire calendar (green dashed lines). 

In many cases, the base simulation captures the general 
shape of the vertical ozone profile quite well (a,b,c,d,f), while 
the simulation with wrong meteorology does not (a,b,d,f)- 
This again illustrates the importance of large-scale transport 
in modeling tropospheric ozone, w-hich we also concluded 
from our comparison to TTOC in Section 6. As expected, 
adjusting the emission parameterization most strongly influ- 
ences the stations over the Atlantic Ocean and African con- 
tinent (c,d,e,f). What is most important to note from the 
comparison however, is that many features of the individ- 
ual sonde profiles are not reproduced by the model. Adjust- 
ing the biomass burning emission parameterization, or any 
of the other parameters we tested, does not appear to lead to 
systematic improvements. 

There are several reasons for this. First of all, we men- 
tion the difference in vertical resolution, which strongly lim- 
its the amount of detail our model can produce. Secondly, 
we mention the strong damping present in the ozone produc- 
tion, which is obvious from the small differences between 
the blue and green lines in Fig. 7 and from Table 3. Gener- 
ally, the ozone mixing ratios increase by approximately 0-5 
ppbv at all altitudes, leading to an increase of several Dob- 
son Units when vertically integrated. This difference is small 
compared to the observed underestimate of TTOC. Finally, 
the comparison presented is obviously hampered by the rep- 
resentation of different spatial scales. When we try to im- 
prove our model through changes in the parameterizations, 
the response seems to be strongest on the large scales which 
are dominated by transport in our model, while the smaller 
scales remain largely unaffected. Since these smaller scales 
are dominant in the ozone sondes, intercomparing the dif- 
ferent model runs does not convey much information on the 
performance of the models parameterizations. In fact, it is 
impossible to judge from Fig. 7 which simulation performs 
better; green or blue? 

The comparison illustrates that individual sondes are not 
the best product to use when testing a CTM. As was recog- 
nized previously by Logan (1999), climatologies of ozone 
observations are more useful in this respect (see also Pe- 
ters etal. (2001)), even though they lack interannual vari- 
ability and the density of the measurement network is quite 
low. We show that the use of TTOC offers a suitable alterna- 
tive, where the lack of vertical information is ’compensated’ 
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Fig. 7. Selected sondes from six SHADOZ stations in 1998. The profiles show observations (black), a standard model run for 1998 (blue 
dashed), a model run with meteorology from 1997 (red dashed), and a model run with additional NO* emissions and a new fire calendar 
(green dashed). 


by the geographical coverage, the frequency of observations, 
and the length of the observational record. 

10 Discussion 

Our model evaluation depends to some extent on the accu- 
racy of the TTOC data as well as their assumed standard 
deviations. The meteorological variability in the model and 
in the satellite TTOC dataset is very similar, which provides 
some confidence in the retrieval algorithm applied to produce 
this dataset. Even though the MR retrieval applies sonde ob- 
servations from two specific years (1990-1991) to produce 
the entire 1979-1992 record, the details of individual years 
seem to have been retained well. Nevertheless, the accuracy 
and the standard deviations of the TTOC observations are not 
yet well quantified because the ozone sondes did not provide 
an absolute reference due to the limited spatial and temporal 
coverage (Thompson and Hudson, 1999). 

To illustrate possible effects, we calculated the values of 
e\ and €2 after applying a uniform negative 5 DU offset to 
the TTOC measurements. This brought model and measure- 
ments almost within one- sigma, and was an even more effec- 
tive way to reduce the underestimate of TTOC than any of 


our sensitivity simulations. However, a 5 DU offset in TTOC 
cannot be justified, and the dependence of 62 on the accu- 
racy of the measured TTOC columns does not invalidate our 
approach. If the TTOC measurements change, for instance 
due to an update in the retrieval algorithm, 62 can easily be 
re-defined to focus on regions that consequently appear most 
problematic, and then applied to the sensitivity runs to isolate 
possible solutions. Nevertheless, a thorough validation of the 
retrieved TTOC data using ozone sondes, as well as a more 
detailed quantification of the uncertainties would add to the 
usefulness of these data in global modeling studies. 

The fact that a 5 DU reduction in the observations would 
have greater impact than any of our sensitivity simulations 
further illustrates the robustness of ozone to perturbations. 
This is also obvious from the small ranges in 61 and in 
Table 2. Although the change in boundary conditions is of- 
ten quite extreme (factor of two), the response of ozone is 
only small. Table 3 illustrates that this nonlinear, damped re- 
sponse is mostly due to strong feedbacks in the NO x -HOx 
chemistry, as well as in the deposition working against any 
change in the concentrations of ozone, NO x , and HNO 3 . 
These feedbacks are especially important in the tropics, since 
both photochemistry and uptake of trace gases by vegeta- 
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tion depend strongly on the sunlight intensity. It is also the 
non-linearity of the system that requires CTMs to use so- 
phisticated mathematical solvers for the chemistry with small 
time steps to solve the system of partial differential equations 
(PDE’s), which generally consumes a large amount of com- 
puter resources. 

An effective way to reduce the demand of computer re- 
sources is by pre-calculating the solution of the system of 
PDE’s for many different boundary conditions and initial 
conditions. The solution of the system at each time step is 
then retrieved from an archive of solutions matching the forc- 
ing. Klonecki and Levy (1997) have shown that this can be 
done in a way that retains many of the feedbacks, and their 
model performs well in simulating tropical ozone, CO and 
NO x . However, the interaction between NO x and HO x was 
not included in their simulations, which removes part of the 
nonlinear response of ozone. Therefore, the method of Klo- 
necki and Levy (1997) produces more ozone over the tropical 
Atlantic and displays greater sensitivity to forcing (see also 
Galanter et al. (2000)). We speculate that the sensitivity of 
such models will usually be larger than in models that calcu- 
late the solution of the PDE*s at each time step. 

Our results show that the meteorology ' particularly 
strongly influences the distribution of ozone. It is striking 
that the range of values of e\ in Fig. 5 is much larger than the 
range in Table 2. Apparently, zonal gradients in TTOC are 
largely determined by transport, and the location of lightning 
NO x emissions that are coupled to transport. Using actual 
meteorology has clear advantages over ’’climatological” me- 
teorology, because the interannual variability is large, also 
addressed by Peters etal. (2001). Although the differences 
in correlation might seem quite low, we have shown that 
they are much greater than can be expected from random 
perturbations. The fact that improvements can be achieved 
indicates that the differences associated with transport are 
an important aspect of the models performance. Neverthe- 
less, many models still rely on climatological wind fields, or 
apply one year of meteorological analyses to simulate any 
year of interest. This approach will yield results that are best 
compared to climatologies of observations, and are much less 
suited to study multi-year time series, or instantaneous dis- 
tributions of species with lifetimes that are smaller than, or 
comparable to, the time scale of transport. 

The requirements of the model concerning meteorology 
are very demanding. Not only does large scale transport de- 
termine the shape of the zonal TTOC distribution, meteo- 
rology also introduces interannual variability, and its inter- 
action with biomass burning emissions is one of the most 
sensitive model parameters. However, meteorological condi- 
tions in the tropics are complex, and their representation in 
numerical weather prediction models such as the ECMWF 
model suffers from a lack of in-situ measurements that are 
used as constraints. Especially subgrid scale vertical mixing 
processes are difficult to model. This includes the vertical 
transport of biomass burning plumes. Mauzerall et al. (1998) 


estimated that the majority of ozone production takes place 
in these plumes. It thus seems obvious that besides improve- 
ments of our knowledge on the occurrence and intensity of 
fires, further improvements should focus on the sub-grid rep- 
resentation of turbulent and convective transport of biomass 
burning plumes. 

Our study shows most convincingly that it is very unlikely 
that a single factor affecting ozone is responsible for the ob- 
served differences between model and observations. Given 
the serious disagreement, and the sensitivities as presented 
in Table 2, a combination of processes would be required 
to bring observations and model within one-sigma. This 
fact, together with the strongly damped response of ozone 
to the applied forcing strongly complicates attempts to in- 
fer information on model parameters (e.g. lightning NO x 
and biomass burning emissions) from the tropical distribu- 
tion of TTOC. The simultaneous use of global CO and NO x 
observations, for instance, could be useful in this respect, as 
these could help to constrain the emissions. Ongoing devel- 
opments in remote sensing might make these products avail- 
able on spatial and temporal scales suitable for the approach 
presented in this work. 

11 Conclusions 

We have combined satellite observations of TTOC with 
model calculations to test the performance of our model and 
prioritize directions for improvement. Such a systematic 
comparison of model calculations with satellite observations 
has not been performed before, and it yields interesting 
insights in problems and possibilities associated with the 
use of increasingly available space-bome tropospheric 
observations. Through a large number of simulations, we 
were able to quantify the influence of meteorology, as well 
as biomass burning and lightning NO x emissions on the 
simulation of TTOC. Our most important conclusions are: 

- With a standard set of emissions such as in Houwel- 
ing etal. (1998); Lelieveld and Dentener (2000); Peters 
etal. (2001) the Atlantic TTOC maximum is not suf- 
ficiently captured, whereas the interannual variability 
over the Pacific Ocean is represented quite well. 

- Using actual meteorology in a model is the most im- 
portant prerequisite in reproducing the horizontal and 
vertical distribution of tropical tropospheric ozone. 

- Efficient feedbacks in the tropical ozone budget strongly 
damp the response of ozone to changes in precursor 
emissions. 

- Systematic improvements in our model can be achieved 
through a better representation of the biomass burning 
timing and magnitude, including their interannual vari- 
ability. Further studies should focus on these processes. 
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as well as the parameterization of lightning NO x emis- 
sions. 

- Quantitative use of satellite derived TTOC to test a 
global model is a promising approach, that would 
strongly profit from a more detailed quantification of 
the uncertainties in TTOC data. 
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